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Abstract 



We will discuss a one-dimensional approximation for the problem of 
wave propagation in networks of thin fibers. The main objective here is to 
i-C . describe the boundary (gluing) conditions at branching points of the lim- 

iting one-dimensional graph. The results will be applied to Mach-Zehnder 
fH , interferometers on chips and to periodic chains of the interferometers. The 

latter allows us to find parameters which guarantee the transparency and 
slowing down of wave packets. 
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1 Introduction. 

The paper concerns the asympototic theory of wave propagation in networks e 

■ of thin fibers when the thickness of fibers e goes to zero. An approximation of 
Q\ | a wave or heat proceses in such a network by a one-dimensional problem on the 

■ limiting metric (quantum) graph has been discussed in physical literature for at 
least three decades. In recent years it was the subject of several mathematical 

• *h . conferences. The central point of the asymptotic theory is the structure of the 

?\ ' physical field near junctions (branching points of the network). In the majority 

of publications on quantum graphs the gluing conditions on the vertices of the 
graph (they correspond to the junctions) have the simplest Kirchhoff's form. 
This form can be justified in some cases (say, for the heat transport in a network 
with insulated walls, pj], [12] ). We discussed some possible applications of 
the quantum graph approximations to the study of periodic optical systems: 
structure of the spectrum, scattering, slowing down of the light, pjI]-[2T]. We 
assumed there that the Kirchoff 's GC at vertices were imposed. 
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However, our recent study leads us to the conclusion that the KirchofPs GC 
are an exception in optical applications where the spectral parameter A is greater 
than the threshold Ao- The latter is equivalent to the condition that the propa- 
gation of waves in the waveguides (cylindrical parts of the network) is possible. 
If fl e is unbounded, this also means that A belongs to the absolutely continuous 
spectrum of the problem. While many particular cases of that problem with 
A = Ao + O(e 2 ) or A < Ao were considered (see [l]-[3Q]), the publications [2"2"]-[25] 
were the first ones dealing with the case A > A , and the first ones where the 
significance of the scattering solutions for asymptotic analysis was established. 
Papers [12]- [5S] contain asymptotic analysis of the spectrum, resolvent and so- 
lutions of the problem in a network f2 £ when e — > 0. It was shown there that 
the GC in those cases have general symplectic structure and can be expressed 
in terms of the scattering matrices defined by individual junctions. 

The main goal of the present paper is to describe the wave propagation 
through networks of thin fibers of necklace type resulting in general GC at the 
vertices of the limiting graph. We say that a network is of necklace type if 
it is periodic in one direction and is bounded in the orthogonal plane. The 
transition from the networks of thin fibers to the one-dimensional problem on 
the graph will be recalled in the next section (see details in 22J-|25J. The 
necklace type graphs will be considered in section 3. We will calculate the 
propagator through one period, find the dispersion relation, describe the band- 
gap structure of the spectrum and find reflection and transmission coefficients 
for the truncated graph. 

The main feature which distinguishes the graph theory from the Bloch theory 
of ID periodic Schrodinger (Hill) operators is that the propagator through one 
period is not an analytic function of the frequency anymore, but a meromorphic 
one. The corresponding poles (resonances) play an important role in applica- 
tions. In the last section, we consider a specific necklace device and show how 
the earlier results allow one to find parameters which provide slowing down of 
the wave packets (slowing down of the light) accompanied by the transparency 
(almost zero reflection) . 

2 Transition from networks to quantum graphs. 

Consider the stationary wave (Helmholtz) equation 

H e u = -Ait = oj 2 u, xeil e , Bu = on 90 E , (1) 

in a domain fl £ c R d , d > 2, with infinitely smooth boundary (for simplicity), 
which has the following structure: Q £ is a union of a finite number of cylinders 
Cj jE (which will be called channels) of lengths lj, 1 < j < N, with diameters 
of cross-sections of order O (e) and domains (which will be called junctions) 
connecting the channels into a network. It is assumed that the junctions have 
diameters of the same order O(e). The boundary condition has the form: B = 1 
(the Dirichlet BC) or B = ^ (the Neumann BC) or B = + a(x), where 
n is the exterior normal and the function a > is real valued and does not 
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depend on the longitudinal (parallel to the axis) coordinate on the boundary of 
the channels. One also can impose one type of BC on the lateral boundary of 
fl e and another BC on the free ends (which are not adjacent to a junction) of 
the channels. 




Figure 1: An example of a domain Q E with four junctions, four unbounded 
channels and four bounded channels. 

The domain Q E shrinks to a one-dimensional metric graph r as e — > 0. The 
axes of the channels form edges Tj, 1 < j < N, of T, and the distances between 
points of Tj are defined by the distances between the corresponding points of 
the channels. The junctions shrink to vertices of the graph T. We denote the 
set of vertices Vj by V. 
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Figure 2: Another example of a domain fi e , network of a necklace type which 
will be studied in more detail later. 

For the sake of simplicity, we impose the following two geometrical condi- 
tions. First, we assume that all the channels Cj i£ have the same cross-section 
7r e (the general case is studied in [53] ). The second condition concerns the 
junctions. We assume that they are self-similar. The latter means that there 
exist an e- independent domain J v and a point x — x v such that 

Jv,e = {{x + ex) : x £ J v }. (2) 

From the self-similarity assumption it follows that ir E is an e-homothety of 
a bounded domain it c R d ~ 1 . 

Let Ao < Ai < A2... be the eigenvalues of the negative Laplacian — A^_i in 
7r with orthonormal eigenfunctions {(fin(y)}, 

-A d -iip n (y) = \„<p n (y), Bu = on dir, 
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where B is the boundary operator on the channels defined in fT}. Then £~ 2 A„ 
are eigenvalues of — A<j_i in tt e and {e~ d / 2 ip n (y /e)} are the corresponding or- 
thonormal eigenfunctions, 

-Ad-xtfniy/e) = s~ 2 X n ip n (y/e), Bu = on <9-7r £ . 

We will call the point £~ 2 Ao the threshold, since it is the bottom of the absolutely 
continuous spectrum of operator ([T]) if il e has an infinite channel. 

We introduce Euclidean coordinates (z, y) in channels Cj_ £ chosen in such a 
way that the z-axis is parallel to the axis of the channel, hyperplane Ry^ 1 is 

Figure 3: Local coordinates in a channel. 

orthogonal to the axis, and C Jje has the following form in the new coordinates: 

Cj. e = {(z,ey) :0<z<lj, yen}. 

If a channel C JiE is bounded (lj < oo), the direction of the z axis can be chosen 
arbitrarily (at least for now). If a channel is unbounded, then z = corresponds 
to its cross-section which is adjacent to the junction. 

We will impose the condition (eui) 2 G (\q, Ai) (see [23] for the general case). 
Note that waves governed by {1} do not propagate through the channels if 
{slu) 2 < Ao- There exists only one propagating mode 

e>*"Mv/e), <? = Vu 2 -e- 2 \ , (3) 
if {suj) 2 G (Ao, Ai), and there are many similar modes 

e ±te 'Vi(w/e), (Tj = ^ 2 -e- 2 \ 3l 

if {euj) 2 > Xj. 

While many particular cases of problem (TT]) with (euj) 2 = Xq + 0(e 2 ) or 
(eui) 2 < Ao were considered (see references in [25]), the publications [22]-[25] 
were the first ones dealing with the case (ecu) 2 > Ao, £ — > 0, and the first ones 
where the significance of the scattering solutions for asymptotic analysis of H £ 
was established. In particular, it was shown there that in both cases (eu>) 2 > Ao 
and [eoj) 2 rs Ao, the scattering solutions for equation (TT]) and the resolvent of 
the operator H £ can be approximated by the corresponding solutions of the one- 
dimensional problem on the limiting graph T with the GC expressed in terms 
of the scattering matrices of the individual extended junctions. 
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Let us recall the definition of scattering solutions for the problem ([TJ) in f2 £ 
when (eui) 2 e (X , Ai). The scattering solution W = ^ p _ £ describes the propaga- 
tion of an incident wave with unit amplitude and frequency uj coming through 
the channel C P)E . 

Definition 1 Let Ao < (e^) 2 < Ai. A function VP = ^ P , 6 , 1 < V < m i * s called 
a solution of the scattering problem in tt E if 

(-A-w 2 )f = 0, ie!! £ ; B* = on dfl E , (4) 

and ^ has the following asymptotic behavior in infinite channels Cj. £ , 1 < j < 
m : 

%,e = [S P , j e- iaz +t pJ e i ' TZ }<p Q (y/e) + 0(e-^), z -> oo, a > 0. (5) 

Here a — y/uj 2 ~ e~ 2 \q, S Pt j is the Kronecker symbol, i.e. S p j = 1 if p = j, 
§ pj =0if p ^j. 

Remark. The term with the coefficient 5 p j in ([5]) corresponds to the inci- 
dent wave (coming through the channel C Pt£ ), t PtP is the reflection coefficient, 
the terms with coefficients t p j,j ^ p, describe the transmitted waves. The 
coefficients t p j — t p j(e,uj) depend on e and uj. The matrix 

T = [tpj] (6) 

is called the scattering matrix. 

Standard arguments based on the Green formula provide the following state- 
ment. 

Theorem 2 When Xq < (ecu) 2 < Ai, the scattering matrix T is unitary and 
symmetric (t Pl j = tj tP ). 

It happens that the scattering solutions ^ PtS can be approximated with an 
exponential in e accuracy using the scattering solutions of a one-dimensional 
problem on the limiting graph T which are defined as follows. Consider the 
following equation on T 

-^V = "V, o^^uj 2 -e- 2 \ . (7) 

Obviously, 

on the edges Tj CT. 

We split the set V of vertices v of the graph into two subsets V = V\ U Va, 
where the vertices from the set V\ have degree I and correspond to the free ends 
of the channels, and the vertices from the set Vi have degree at least two and 
correspond to junctions. 
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Definition 3 We will say that %j) — tp p ^ e is a solution of the scattering problem 
on the graph T with the incident wave coming through the edge T p ifipp,e satisfies 
equation 

VpA-y) = S P^ az + tp,j£ iaz , 7 G r„ (8) 

on infinite edges < j < m, and satisfies the following GC at vertices v of 
T: 

Bip = atveVi, (9) 
i[I v + T v {eu)}^ v \z)-a[I v -T v {eu)}^ v \z) = Q, z = 0, v € V 2 . (10) 

We keep the same BC at v G V% as at the free end of the corresponding channel 
of fl e , see ([1]), and we will specify GC (|10[) in the next two paragraphs. However, 
first we would like to stress that scattering coefficients t p j in (|8|) are not required 
to coincide with those defined in ((5]) for the problem in the domain fl £ . With 
the appropriate choice of GC (|10|). these coefficients are the same, and therefore 
we use the same notation. 

We choose the parametrization on T in such a way that z = at v for all 
edges adjacent to this particular vertex. Let d — d(v) > 2 be the order (the 
number of adjacent edges) of the vertex »eF 2 - For any function tp on T, we form 
a column- vector ifj^ = ijj^ v \z) with d(v) components which is formed by the 
restrictions of -0 on the edges of T adjacent to v. We will need this vector only for 
small values of z > 0. The GC (|10[) are defined in terms of auxiliary scattering 
problems for extended junctions J^° e . Each extended junction consists of J„ jE 
and all the channels adjacent to J v , e . If some of these channels have finite 
length, we extend them to infinity (see Fig. [4|). The matrix T = T v (euj) is 



Figure 4: An extended junction for a necklace waveguides. 

the scattering matrix for the problem (fTJ) in J^° £ and I v is the unit matrix of 
the same size as the size of T. Note that the self-similarity of J^° e implies that 
T = T v (ew) depends only on the product eu. Hence, T v can be determined 
by solving the scattering problem in the corresponding extended junction with 
e = 1. Now (ITU1) is defined. We need only to take components of the vector ipw 
in the same order as the order of channels of . 

Denote by F(e) the set {<jJj} of the values of u> for which problem ([7]), ©, 
(TTU)) has a nontrivial solution ip £ L 2 (T) and (sujj) 2 G [Ao,Ai]. Note that both 
equation and the boundary conditions of the problem depend on to and e. Let 
F° be the set {ujj} of values of u such that uj? is an eigenvalue of the operator H e 
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(see (P)) in one of the domains J£° £ and (euij) 2 € [A , Ai]. Due to self-similarity 
of the extended junctions J£° £ , the set eF° does not depend on e. Let F v be 

the e~-neig hborhood of the set F° U F(e). 

Theorem 4 (\2S§) 1) The set F° U F(e) has finitely many points (the number 
of points depends on e). 

2) For any interval [Ao, A'], A' < Ai, there exist p, v > such that scattering 
solutions typ^ e (x) of the problem in Q e have the following asymptotic behavior 
on the channels of fl E as e — ► 

%, e (x) = ^ e) (7)vo(|) +r { p £) (x), 7 G r, 
where ipp (7) are the scattering solutions of the problem on the graph T and 

\rV(x)\ < Ce^ 1 , M 2 e [A ,Y], w i F\ 

Here 7 = y(x) is the point on T which is defined by the cross-section of the 
channel through the point x, and d(^) is the distance between 7 and the closest 
vertex of the graph. 

Note that this theorem implies the coincidence of the scattering matrices of 
the problems on f2 e and on the graph T. 

We will conclude this section by the following important proposition: 

Proposition 5 Suppose that det[I v + T v (eui)] is not identically equal to zero. 
Then for all u> such that (sto) 2 G [Ao, Ai], except at most a finite number of 
points, the GC blCfy can be written in the form 

j t ^\t)-aA v ^ v \t)=Q, t = 0, o 2 = ^uj 2 - £- 2 A , v G V 2 , (11) 
where matrix 

A v = -i[I v + T^eu)]- 1 ^ - T v {eiu)} 
is real valued and symmetric (A v = A' v ). 

Proof. Due to analyticity of T v in A, we need only to justify properties of A v . 
From Theorem [2] it follows that the eigenvectors of T v can be chosen to be real 
valued, i.e., there exists a real- valued orthogonal matrix C v — C V (X) such that 
T v = C~ 1 D V C V where D v is a diagonal matrix with elements Zj, \zj\ = 1 on 
the diagonal. Numbers Zj are (complex) eigenvalues of T v . Obviously, 

I v iT v = I v dzC v D V C V — C v [I v iD v ]C V: A v = iC v [I v -\-D v ] [I v D V \C V . 

It remains to note that the numbers i ], Zj , are real. ■ 
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3 Necklace graphs 



In this section we will consider a periodic metric graph T of the necklace type 
(see [20]). It has the following form: one cell of periodicity consists of two arches 
Ci,C% of lengths l\ and I2 < h connected at end points and of a segment £3 
of length Z3 starting at one of these points (see fig. 1). We assume that the 
necklace is placed horizontally. Let ^m-b^m be the left and the right end 
points of the segments with vq being the origin, and with two arches connecting 
the points i>2m and W2m+i, m = 0, ±1,±2, ... . 

We introduce two related local coordinates z and s on the edges of the graph. 
Both are the lengths of the corresponding part of an arch or a segment, measured 
from some end of the edge. When a neighborhood of some vertex is considered 
(for example when GC are defined), the distance is measured from that vertex 
for all the edges adjacent to this vertex. The coordinate (distance) z is used in 
this case. In other cases it will be convenient for us to measure the distance 
from the left end of the edge to the right. We will specify this situation by using 
parameter s instead of z. Thus, s = z or s = I — z where / is the length of the 
edge. 

We equip the graph with the natural Lebesgue measure and consider the 
Hamiltonian H on L 2 (T) given by H = —^2 = ~j^2 on the graph with the 
following GC (see (HU) 

j t ^ v) (z) - oA v (ecj)^ v) {z)=0, z = 0, v&V, (12) 

at the vertices of T. Here and below we use 

a = V 'lu 2 — e~ 2 \q. 

Our goal in this section is to study the propagation of waves on T governed 
by the equation 

Hi' = -^ = a 2 ^, 7 er, (13) 

and GC (H2). Note that both the equation and GC depend on the frequency oj 
and e. 

We enumerate the components ipj of the vector ip( v ) in the following order: 
ipj corresponds to the edge of the length lj . Thus the first two components of 
the vector correspond to the shoulders of the loop, and the third component 
corresponds to the straight edge. 




Figure 5: Graph T for the necklace waveguide. 
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The goal of this section is to define and evaluate the (Priiffer) monodromy 
operator M a (transfer operator over the period) for the problem (fT3|) , (| 12[) . Let 

ip' = £^ and ( -1.0' 1 ( a ) ^ e ^ c Cauchy data (it always will have the factor 

cr -1 in the second component) of the solution i/j of the equation (fT3)) evaluated 
at the point a of a straight segment of T. When a — v n is one of the end points 
of a segment, we understand this vector as the limit of the corresponding vectors 
evaluated at a as a approaches v n moving along the segment (not along one of 
the arches). We denote by M a the monodromy operator: 

M„ : ( a % ) M - ( a % ) M 
and we denote by T a the (Priiffer) transfer operator over the loop: 

a"V ) M ^ ( ) (Vl) - 

We will use the same notations M a , T a for the matrices of the operators as for 
the corresponding operators. 

Let us write matrix A v = (aij) (see (fTT|) . (fT2|l) in the form 

)•«-(* HE)—* 

(14) 

We will need the following matrices 

* inah ?),C=(™ ah ,° V P = C + SB, (15) 
sin erf 2 J \ cos al 2 J 

M = (I-P 2 )- 1 PS, N = (I — P 2 )~ 1 S. (16) 



Theorem 6 1) Matrix T a has the form 

/ml \ 

T* = -( A»" ro ), n = (S,NS), m = c+(S,M5). (17) 
2) Matrix M a has the form 

( cosa/3 sina/s \ 
\ — sm«3 cosct( 3 y 

Proof. Let i/'j = tyj^fi), j — 1,2, be the restrictions of the function i\> on the 
upper and lower arches of the graph between the points vq and v\. Obviously, 

, , ^(0) sincr(Z 1 — s) +tpi(li)smas 

Vi{s) = : j • 

smCTii 
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Thus, 



-o-ipi(0)cosah+(xipi(h) ,, -atpi(0) + aipi(h) cos ah 
Vi(0) = -r—, , V'xC'l 



sin cr^i sinali 

A similar formula is valid for ip2 ■ Hence, taking into account the relation between 
parameters s and t, we obtain the following connection between the values of 

the vector <fi — f ^ ~\ and its derivative ^ at points vq and v\ of the graph: 

^(v Q ) = -aCS-^M + aS' 1 ^), = -aCS- 1 ^) + aS' 1 ^). 

dz dz 

Let us denote the restrictions of ip on the straight edges of the graph by ipz ■ 

Consider vectors V ( " o) = ( f^°\ ) , ^ (vi) = ( f ^\ ) . GC STR implies 

-C5-V(«b)+5-V(»i) AM a \ / 0(«o) 

-C5-V(»i)+5-V(«b) \ — { B S\( 4>{vi) 
tr-^Mvi) ) V <** c / V ^a("i) 

These equations can be rewritten in the following form 

-(CS- 1 + B)cj>{v ) + STV(«x) = #3(«o), 

a^^-^(vx) = 8*4>(vi) + cip 3 (vi). 
as 

We multiply the first and third equations by S, replace C + SB by P (see (fTS)) 1 ) 
and then solve these equations for (j)( v o), 4>{ v i)- This implies 

= -(/ - P i> )- 1 PS6tI>s(vo) + PY'SSMvi), 
<f>(vi) = -(I - P 2 )~ 1 S5i/j 3 (vo) + (I - P 2 )" 1 ^^^). 

We substitute these relations into the second and forth equations of the system 
above and then solve those equations for ip3(vi), ■4-ip3(vi). This provides 
the transfer operator over the loop with the transfer matrix defined in (|17p . 
This completes the proof of the first part of theorem. The second statement of 
the theorem is obvious, since the left factor in the right-hand side of (fT5| is the 
transfer matrix over the segment [t?o,Ui] of the graph. ■ 



10 




4 Some applications of necklace waveguides 

We will discuss here two practical features of necklace waveguides. The first con- 
cerns slowing down of the light ( slowing down of propagation of wave packets) 
in these waveguides. There is an extended literature on the principles of this 
phenomenon, possible applications and practical devices. Usually some type of 
a periodic structure is suggested for these devices with a band-gap structure of 
the spectrum. If a narrow band is created, then the dispersion relation for the 
corresponding frequency is flat and the group velocity V g is small. 

The main feature which distinguishes a necklace waveguide from other one- 
dimensional (or quasi one-dimensional) periodic problem is the following. While 
the Hill discriminant for a standard periodic Schrodinger operator is an analytic 
function of frequency, it is meromorphic for the necklace waveguides. We will 
show that one can easily find parameters when the Hill discriminant has two 
close poles (as close as one pleases) with a zero in between. Thus the band 
will be as narrow as we please around a chosen value a = Co of the frequency. 
Therefore, the group velocity will be small if the support of the wave packets 
belongs to a small neighborhood of ctq. 

The second feature concerns the truncated necklace graph (waveguide) Tn 
which consists of N cells of periodicity confined between points vq and V2N and 
the rays (—00, vq), (v2n,oo.) When propagation of a narrow in frequency wave 
packet through a finite device is considered, it could happen that the waves, 
which are slowing down in the device, can not enter the device. In other words, 
one needs to know that the reflection coefficient r for the truncated necklace 
waveguide is not too big (better if it is small) for the frequencies where the 
slowing down occurs. It will be shown that for a given frequency a = ao one 
can choose the parameters lj, j — 1, 2, 3, such that V g = r = at a = ao, and 
therefore they are small for a close enough to <tq. 

Slowing down, preliminary discussion. Theorem 0] allows us to reduce the 
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study of propagation of single frequency waves and wave packets in periodic 
necklace waveguide fl e (see Fig. [2} to a study of the corresponding problem on 
the necklace graph T (see Fig. [5]). One needs only to take in correspondence 
the frequencies ui and a of the waves in Q e and on T using the relation 



(7 = \J\jJ 2 - £" 2 A . 



The spectrum of periodic problem (|13p . (|12p has a band-gap structure with 
the bands on the er-axis defined by the inequality 

\F(a)\ < 2, F(a) = TrM a , 

where function F (called the Hill discriminant) is equal to the trace of the mon- 
odromy matrix (see Fig [5]). The same function defines the dispersion relation 
(we will write it in the form k = k(a)) of the problem on T: 

cosk(a) = ^F(a) = ^TrM a . (19) 

Let lu = u>o belong to the frequency interval of a narrow wave packet of the 
problem in Q £ and 



cr = y - £ 2 V 

We will find parameters lj in such a way that 

in^o)| = 0, (20) 

and there is a point a = a\, for example, to the right of <To, such that F{a) has 
a pole at tr = <j\ and |<7o — a%\ = e <S 1. Then there is a point a' G (cfq, a\) such 
that 

|F(a)|<l for ae A = (a ,a'), \F(a')\ = 1. 

Then interval A belongs to a band. From (JTHJ) it follows that k(a) changes by 
7r/6 between <jq and a'. Hence, k'(a) = 0(l/e) on A or on some part A' of A, 
i.e., 

L 
k~a) 

Note that matrix T a has a pole if n — 0. Since matrices T a and M a differ 
by a rotation, one can easily choose I3 such that F(a) = ^TrM CT has a pole 
at a point where n = (in fact, the pole of F(a) exists for all but specific 
values of Z3). Thus the following equation provides the poles of F(a) under an 
appropriate choice of l 3 (see (fI5 ]) -([T7 )) ) 

n = (5, (I - P 2 )- 1 ^)) = )-{8, [{I - P)- 1 + (7+ P^SS)) =0, P = C + SB. 

Let us introduce the matrix 

T= I —^=( X n °). x = Un^, y = tan^. (22) 
S \ y )' 2 ' J 2 v 1 



V 9 = — = ^^ ff eA'. (21) 
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Then the equation n = can be rewritten in the form 

2n = (8, [(T - B)- 1 + (T- 1 + B)^ 1 }8)) = 0. (23) 

We will conclude this subsection by an re-writing (|23|) using the components of 
matrix B and vector 8 (see (fT4]0 : 

r, (y~ 0.2,2)81 + (x - ai i)<5| + 2ai i2 <5i52 
(x - ai,i)(2/ - 02,2) - a\ 2 

, (y _1 + "2,2)^ + (x~ x + a M )<5f - 2a li2 ^^ 



(x 1 +ai,i)(y 1 + 02,2) - a\ ■ 



= 0. (24) 



This is an algebraic equation of the forth order with respect to (x, y). We will 
discuss it in more detail later. 

Transparency, preliminary discussion. Let us recall the estimate (see |21| ) 
for the reflection coefficient r — r^ by a finite slab of periodic media which 
consists of N periods: 

sm^) 

where H-Mojl is the Gilbert-Shmidt norm of the monodromy matrix, i.e. 



||M a H 2 = ||(mij)l| 2 = £ 

i,J<2 



We will choose parameters lj in such a way that the frequency support A of the 
wave packet is in the middle of a band where ^(c)! < 1. Then | sinfc(cr)| > ^ 
there, and 

\r N (a)\ = A|(||A/ (T || 2 -2) 1 / 2 , a E A. (26) 

We will choose h, I2 in such a way that ||M CT[) || 2 = 2 , and therefore 

r N (a ) = 0. (27) 

Then rAr(cr) is small in a small neighborhood of o~o 7 and we will have both the 
slowing down of the light and the transparency in A'. 

Note that \\M a \\ = ||T CT || (see JIBJ)) and that detT,, = 1. Thus ||T CT || > 2 
and 1 1 1 1 = 2 if matrix T a is orthogonal (the sum of its off-diagonal elements 
is zero). Hence, condition 

(to + n)(m - n) + 1 = to 2 - n 2 + 1 = (28) 

provides (|2"T|) . We are going to write this condition more explicitly. 
From fIS]). (HO it follows that 

m+n = c+{8, (I-P^SS) = c+{5, [S'^I-C-SB)]-^) = c+{8, [(T-B)]-^), 
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where T is defined in (J22). Similarly 



m- n = C -(S, (I+py^S) = c-(S, [S^il+C+SB)}-^) = c-(6, [{T^+B)]-^). 

We substitute the last two formulas into (f2"5| and obtain the following exact 
form for 

r . (y - a 2 ,2)$l + (x - ais)Sj + 2ai i2 Si$2 , 



{x - ai.i)(y - a 2 .2) - a\a 
(?y _1 + a 2 . 2 )S 2 1 + (x- 1 + oi,i)tf| - 2a h2 S 1 S 2 , 



(x 1 + ai,i)(y 1 + a 2>2 ) 



1 = 0. (29) 



Thus, the set of transparency points {rpj — 0} is also given by zeroes of a 
polynomial of forth order in (cc, y)-plane, x = tan y = tan To be more 
accurate, we need to omit points where n = from this set, since (|2"51) provides 
the orthogonality of T a only if n ^ 0. 

The choice of parameters. The center ao hi the frequency interval of the 
wave packet is given. We need to choose h, l 2 in such a way that (f27|) -(j2"9 j) 
hold at Co and this point is close to a point where (|24|) . (|23|) hold. Note that 
condition (|27[) is equivalent to the orthogonality of matrix T ao . Then one can 
easily transfer it by rotation to a matrix with zero trace, i.e., one can find I3 
such that (|20"f holds. Hence, it remains to choose l\, h appropriately. 

We note that equation (f2"5| for cr contradicts the condition n = imposed 
by (f2"3"| . Moreover, n can not vanish at a point close to cto if n is smooth. 
In order for equation (|28[) to be valid at a point ctq and n to be zero at ui which 
is close to Co, function n has to be singular near o-q. Thus one must choose 
parameters (x,y) near the point where one of the denominators in (|24[) is zero. 
Let us choose the denominator of the first fraction. Then the numerator of the 
first fraction also must be small (otherwise (|24|) is not valid at a = ci). By 
equating both the numerator and the denominator of the first fraction in (124ft 
to zero we find that 

S 2 

x « a M - ai,2i-, y « a 2 .2 - a i,2-r- 
b 2 b\ 

Finally, we fix small e and choose 

x — ai 1 - ai, 2t~ + e- (30) 

Then we solve (|2T)1) for y in a neighborhood of the point 02,2 — 0,1,2^-- This 
implies 

V = a 2 , 2 - ajj- - (^) 2 e + 7 e 2 + 0(e) 3 , (31) 
01 01 

where the exact value of 7 can be easily found from (|29|) . We put a = uo in 
([22]) and determine Z 2 from (|30|). ([3T1). This choice of J x , Z 2 implies ((27)) . 
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By solving (f2"4"| asymptotically, we find that its solution also has the form 
(j30|) , (f3"Tj) with a different value of 7. This justifies the existence of a pole of the 
Hill discriminant at the distance 0(e 2 ) from the point <tq. In fact, let us omit 
the cubic term in (|3ip. Then equations (fBT))) , (f5T|) with two different values of 
7 define two parabolas Pi and Pi with the same vertex and tangent line at the 
vertex. If we put a — <7q in (|22|) and fix e, we get a point on Pi at the distance 
of order e from the vertex which defines h,fo- This choice implies ((27)) . When 
a changes: a = <7q + t, < t < to, point (a;, y) moves in the direction of the 
vector (I1J2) or in the opposite direction (see[22]). We may change the sign of r, 
if needed, to guarantee that the point moves toward the second parabola where 
n = 0. The only unacceptable situation is when the vector (li, I2) is tangent to 
Pi. It will not happen in a generic case. Besides one can always avoid it by 
changing li or I2, since they are defined up to an integer multiple of 2n. 

The arguments above prove that wave packet with frequencies in 0(e 2 ^semi- 
neighborhood of do will propagate with the group velocity V g =0(e 2 ) and the 
reflection coefficient will have order 0(e 2 ). 

We conclude this subsection by the following remark. Let the geometry of 
the junctions be chosen. After that, we must choose specific I1J2 satisfying 
(1501) . (EH)- ^ i s eas y to choose one of these parameters (for example, ^i) by 
changing the distance between two neighboring junctions, but I2 will be defined 
after that by the geometry of the network. In fact it is not the geometrical, but 
only the optical length which plays role here. One can preserve the geometry of 
the network and change the refraction index in some of the channels to satisfy 
the relations ([30]) . (f3"Tj) . This produces the same effect as changing the lengths 
of the corresponding channels. 
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